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We give a brief (and highly incomplete) overview of the current experimental and theoretical status of high 
energy cosmic rays and their secondary 7— rays and neutrinos. We focus on the role of large scale magnetic fields 
and on multi-messenger aspects linking these three channels. We also recall that the flavor composition of neutrino 
fluxes from astrophysical sources contains information on both the source conditions and neutrino physics. 



1. Introduction 

After almost 90 years of research, the origin of 
cosmic rays is still an open question, with a degree 
of uncertainty increasing with energy [1]. Above 
~ 100 MeV the CR spectrum exhibits little struc- 
ture and is approximated by broken power laws 
oc iT 7 : At the energy E-4x 10 15 eV called the 
"knee" , the flux of particles per area, time, solid 
angle, and energy steepens from a power law in- 
dex 7 ~ 2.7 to one of index ~ 3.0. The bulk of the 
CRs up to at least that energy is believed to orig- 
inate within the Milky Way Galaxy, typically by 
shock acceleration in supernova remnants. These 
objects have been seen in 7—rays up to ~ 10 TeV 
by experiments such as HESS [f] and MAGIC [3], 
with a spectrum roughly scaling as E~ 22 . This is 
consistent with 7— ray production in interactions 
with the ambient gas of primary cosmic rays with 
a similar spectrum at the source. The difference 
to the cosmic ray spectrum observed to scale as 
E~ 2J below the knee can be explained by diffu- 
sion in the galactic magnetic fields. 

Above the knee the spectrum continues with a 
further steepening to 7 ~ 3.3 at E ~ 4 x 10 17 eV, 
sometimes called the "second knee". There are 
experimental indications that the chemical com- 
position changes from light, mostly protons, at 
the knee to domination by iron and even heav- 
ier nuclei at the second knee [3). This is in fact 
expected in any scenario where acceleration and 
propagation is due to magnetic fields whose ef- 



fects only depend on rigidity, the ratio of charge 
to rest mass, Z/A. This is true as long as en- 
ergy losses and interaction effects, which in gen- 
eral depend on Z and A separately, are small, as 
is the case in the Galaxy, in contrast to extra- 
galactic cosmic ray propagation at ultra-high en- 
ergy. Above the so called "ankle" or "dip" at 
E ~ 5 x 10 18 eV, the spectrum flattens again to 
a power law of index 7 ~ 2.8. This latter fea- 
ture is often interpreted as a cross over from a 
Galactic component, which steepens because cos- 
mic rays are not confined by the galactic magnetic 
field any more or because Galactic sources do not 
accelerate beyond the ankle, to a harder compo- 
nent of extragalactic origin. However, the dip at 
i?~5x 10 18 eV could also be explained by pair 
production by extra-galactic protons, if the extra- 
galactic component already starts to dominate 
below the ankle, for example, around the second- 
knee [5] at a few times 10 17 eV. This requires a 
relatively steep injection spectrum oc £?- 2 - 6 - 2 - 7 . 
Below a few times 10 17 eV this extra-galactic 
component would become unobservable at Earth 
due to diffusion in extra-galactic magnetic fields 
(EGMF) [6 . In addition, the effective volume- 
averaged injection spectrum has to become flatter 
somewhere below ~ 10 18 eV in order to avoid ex- 
cessive power going into cosmic rays and to avoid 
overproduction of GeV-TeV 7—rays from pp in- 
teractions with the ambient gas. 

The low cross-over scenario also requires the 
dominance of protons around the dip. Theoret- 



2 



ically, this can be achieved either because pref- 
erentially protons are accelerated or because ex- 
tended EGMF lead to strong photo-spallation 
during propagation [7]. Experimentally, above 
~ 10 17 eV the chemical composition is basically 
unknown 8J. Around 10 18 eV the situation is 
particularly inconclusive as HiRes [9] and HiRes- 
MIA [I0 a data suggest a light (proton dominated) 
composition, whereas other experiments indicate 
a heavy composition 4 . In any case, the cos- 
mic ray flux should be extra-galactic at least 
above the ankle, where a galactic origin would 
predict an anisotropy toward the galactic plane 
because galactic magnetic fields can no longer 
isotropize the cosmic rays. No such anisotropy 
is seen. There are also experimental indications 
for a chemical composition becoming again lighter 
above the ankle, although a significant heavy 
component is not excluded and the inferred chem- 
ical composition above ~ 10 18 eV is sensitive to 
the model of air shower interactions and con- 
sequently uncertain presently [8]. In addition, 
should a substantial heavy composition be ex- 
perimentally observed up to the highest ener- 
gies, some sources would have to be surprisingly 
nearby, within a few Mpc, otherwise only low 
mass spallation products would survive propaga- 
tion [llj . In the following we will restrict our dis- 
cussion on extra-galactic ultra-high energy cosmic 
rays (UHECRs). 

No conclusive picture for the nature and dis- 
tribution of the sources emerges yet naturally 
from the data [12 : Arrival directions are approx- 
imately isotropic [15] . suggesting a large num- 
ber of weak or distant sources. But there are 
also indications which point more towards a small 
number of local and therefore bright sources, 
especially at the highest energies: First, the 
AGASA ground array claimed statistically signif- 
icant multi-plets of events from the same direc- 
tions within a few degrees [14113] , although this is 
controversial [T5] and has not been seen so far by 
other experiments such as the fluorescence exper- 
iment HiRes [16] . The spectrum of this clustered 
component is oc E~ 1,s and thus much harder than 
the total spectrum [14]. Second, nucleons above 
~ 70 EeV suffer heavy energy losses due to photo- 
pion production on the cosmic microwave back- 



ground (CMB) — the Greisen-Zatsepin-Kuzmin 
(GZK) effect [17] — which limits the distance to 
possible sources to less than ~ 100 Mpc [18] . This 
predicts a "GZK cutoff" , a drop in the spectrum, 
whose strength depends on the source distribu- 
tion and may even depend on the part of the sky 
one is looking at: The "cutoff" could be mit- 
igated in the northern hemisphere where more 
nearby accelerators related to the local superclus- 
ter can be expected. Apart from the SUGAR 
array which was active from 1968 until 1979 in 
Australia, all UHECR detectors completed up to 
the present were situated in the northern hemi- 
sphere. Nevertheless the situation is unclear even 
there: Whereas a "cut-off" is consistent with the 
few events above 10 20 eV recorded by the fluores- 
cence detector HiRes [19] and with the first data 
release of the Pierre Auger observatory [20] , there 
is a tension with the 11 events above 10 20 eV de- 
tected by the AGASA ground array [21]. Still, 
this could be a combination of statistical and sys- 
tematic effects [22], especially given the recent 
downward revision of the energy normalization 
in AGASA [23] . The solution of this problem will 
have to await more analysis and more statistics 
and, in particular, the completion of the Pierre 
Auger project [20] which combines the two com- 
plementary detection techniques adopted by the 
aforementioned experiments and whose southern 
site is currently in construction in Argentina. Fi- 
nally, about 1% of the HiRes steroa events around 
10 19 eV seem to correlate with active galaxies of 
the BL Lac on a scale of ~ 0.6°, with a signif- 
icance of ~ 10 -4 [24125] . Due to deflection of 
UHECR in the galactic magnetic field this, how- 
ever, would have to be neutral primaries that can- 
not be created in the necessary quantities over the 
distances involved. 

2. Role of large scale magnetic fields 

The hunt for UHECR sources is further com- 
plicated by the presence of large scale cosmic 
magnetic fields which may significantly deflect 
charged cosmic rays even at the highest ener- 
gies, in particular if sources correlate with high 
magnetic field regions such as galaxy clusters. 
A major issue in UHECR propagation studies 
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is, therefore, the strength and distribution of 
EGMF. Galaxy clusters harbor magnetic fields of 
fj,G strength, but it is poorly known how quickly 
these fields fall off with increasing distance from 
the cluster center. The current data indicate that 
fiG strength magnetic fields extend to at least ~ 1 
Mpc [26] and possibly to larger distances |27|28j . 
Beyond ~ 1 Mpc from a cluster core, however, 
probing the magnetic fields becomes extremely 
difficult because the Faraday Rotation Measure 
loses sensitivity in low density regions. Further- 
more, the intracluster magnetic field topology is 
also poorly known, although the situation will 
likely improve in the future, for example with 
the advent of powerful radio astronomical instru- 
ments such as the square kilometer array. 

One possibility in the meantime is to adopt 
large scale structure simulations (LSS) which in- 
clude magnetic fields. In Ref. [29] , the authors use 
magnetic fields derived from a cosmological LSS 
with magnetic fields generated at the shocks that 
form during LSS formation, whereas in Ref. [30j 
and Ref. [31] fields of "primordial" origin have 
been considered. While the different models for 
initial magnetic seed fields produce different large 
scale magnetic field distributions and, therefore, 
lead to different predictions for UHECR deflec- 
tion, there is still a significant discrepancy be- 
tween Ref. [29|30j and Ref. ;3jj. hinting that 
other technical reasons may play a role here. In 
the more extended fields from the simulations of 
Refs. [29130] deflection of protons up to 10 20 eV 
can be up to tens of degrees, whereas deflections 
in the simulations of Ref. [31] are typically below 
a degree. Assuming the EGMF correlates with 
the infrared luminosity density, Ref. [35] recently 
found results closer to Refs. [29130] than Ref. [3"T] . 

We recall that since acceleration is rigidity de- 
pendent, at the acceleration sites the highest en- 
ergy cosmic ray flux is likely dominated by heavy 
nuclei. If this is indeed the case, it is interest- 
ing to point out that even in the EGMF scenario 
of Ref. [3Tj . deflections could be considerable. In 
contrast to the contribution of our Galaxy to de- 
flection which can be of comparable size but may 
be corrected for within sufficiently detailed mod- 
els of the galactic field, the extra-galactic con- 
tribution would be stochastic. Statistical meth- 



ods are therefore likely to be necessary to learn 
about UHECR source distributions and charac- 
teristics as well as EGMF. For example, a sup- 
pressed UHECR arrival direction auto-correlation 
function at degree scales, rather than pointing to 
a high source density, could be a signature of ex- 
tended EGMF [29]. 

Finally, EGMF can considerably increase 
the path-length of UHECR propagation and 
thus modify spectra, especially from individual 
sources, as well as the chemical composition ob- 
served at Earth [33] . 

3. Multi-messenger approach: Secondary 
gamma-rays and neutrinos and their fla- 
vor composition 

The physics and astrophysics of UHECRs are 
also intimately linked with the emerging field 
of neutrino astronomy [34j as well as with the 
already well established field of 7— ray astron- 
omy [35] . Indeed, all scenarios of cosmic ray ori- 
gin, from the galactic scale [36] to top-down [37] 
and Z-burst models at the highest energies |40 a , 
are severely constrained by neutrino and 7— ray 
observations and limits. This has, for exam- 
ple, important consequences for theoretical pre- 
dictions of diffuse fluxes of extragalactic neutri- 
nos above about a TeV whose detection is a ma- 
jor goal of next-generation neutrino telescopes: 
If these neutrinos are produced as secondaries 
of protons accelerated in astrophysical sources 
and if these protons leave the sources and con- 
tribute to the UHECR flux observed, then the 
energy content in the neutrino flux can not be 
higher than the one in UHECRs, leading to the 
so called Waxman-Bahcall bound for transparent 
sources with soft acceleration spectra [41142] , If 
one of these assumptions does not apply, such 
as for acceleration sources with injection spectra 
harder than E~ 2 and/or opaque to nucleons, or if 
much fewer nucleons than 7— rays and neutrinos 
are produced, such as in top-down scenarios, the 
Waxman-Bahcall bound does not apply, but the 
neutrino flux is still constrained by the observed 
diffuse 7— ray flux in the GeV range. 

Fig. [1] provides a sketch of "realistic" cosmic 
ray, 7— ray, and neutrino flux predictions in com- 
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Figure 1. Model fluxes compared to experimen- 
tal data, limits and sensitivities. Primary cos- 
mic ray fluxes (data and a model, see text) are 
shown in black, the secondary 7— ray flux ex- 
pected from proton interactions with the CMB 
and infrared background in red and the "guaran- 
teed" neutrino fluxes per neutrino flavor in blue: 
atmospheric neutrinos, galactic neutrinos result- 
ing from cosmic ray interactions with matter in 
our Galaxy [33], and "GZK" neutrinos resulting 
from cosmic ray interaction with the CMB and 
infrared background. The GZK neutrino fluxes 
depend on the distribution of the (unknown) pri- 
mary cosmic ray sources for which we assumed 
active galactic nuclei (AGNs) above f 17 eV, us- 
ing our public CRPropa code [38139] . Cosmic ray 
interactions within these sources can also pro- 
duce neutrinos for which one example is given 
(AGN v) [40 . The flux of atmospheric neutrinos 
has been measured by underground detectors and 
AMANDA. The dashed and dotted blue lines are 
existing upper limits and future sensitivities to 
diffuse neutrino fluxes from various experiments, 
respectively [34], assuming the Standard Model 
neutrino-nucleon cross section extrapolated to the 
relevant energies. The maximum possible neu- 
trino flux would be given by horizontally extrap- 
olating the diffuse 7— ray background observed by 
EGRET [44]. 



parison with experimental observations, limits, 
and sensitivities. It shows a theoretical scenario 



in which extra-galactic cosmic ray sources roughly 
evolving as quasars inject a spectrum oc E~ 26 
of dominantly protons down to ~ 1 17 eV where 
a cross-over to galactic cosmic rays occurs [5]. 
The " cosmogenic" neutrino flux produced by pro- 
tons interacting with the low energy photon back- 
ground considerably depends on these assump- 
tions which can thus be used to test them [45] . 

Apart from cosmogenic neutrinos produced 
during propagation of UHECR, neutrinos can 
also be produced within astrophysical sources 
such as AGNs (see Fig. [T]) or 7— ray bursts. In 
the absence of matter effects, a source at cos- 
mological distances injecting neutrino fluxes with 
a flavor ratio oc wp, [3 = e, fi, r, leads to a fla- 
vor mixture <f> a oc Yp i w p\U a i\ 2 \Upi\ 2 observed 
at Earth, where U a i is the mixing matrix and i 
labels mass eigenstates. Therefore, if both pions 
and muons decay before loosing energy around 
the source, w e : : w T ~ I : 2 : and thus 
4> e : (fift : <fi T ~ I : I : I. At high energies the 
meson and muon energy loss time t\ oss (E) be- 
comes shorter than their decay time Er/m, and 
the neutrino spectrum will be suppressed by a fac- 
tor ~ mt\ oss (E) I (tE) compared to primary inter- 
action rates. For hadronic cooling, ii oss ^const, 
whereas for radiative cooling at the highest ener- 
gies, t\ QSS (E) oc -E 1 , resulting in a steepening of 
the neutrino spectrum by a factor E~ x and E~ 2 , 
respectively [46] . In addition, at a given energy, 
charged pions decay about hundred times faster 
than muons. There can thus be an energy range 
at which pions but not muons decay before loos- 
ing energy such that w e : u> M : w T ~ : I : 
and thus <j) e : M : 4> T ~ I : 2 : 2. Also, pp inter- 
actions produce both pions of both charges and 
thus give a higher fraction of i7 e compared to pry 
interactions. The observed flavor ratios can thus 
depend on energy and carry information on the 
source conditions [47], but also about the mixing 
matrix itself [4"8] . 

Finally, flavor ratios can probe new physics, 
such as neutrino decay and quantum decoher- 
ence 49 : If all but the lightest mass eigenstate 
j decay before reaching the observer, the flux of 
flavor a observed at Earth would be oc | U a j | 2 , in- 
dependent of the flavor ratio at the source. For 
j = I (normal mass hierarchy) this gives <j) e : 
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0/* : 0-r — 6:1:1, whereas for j = 3 (inverted 
mass hierarchy) one has (/> e : 0^ : T ~ : 1 : 1, 
which should be easy to distinguish from the nor- 
mal case. This would allow to probe lifetimes 
of the order r/m ~ 300(E/ PeV)- 1 a/eV, which 
could improve on current limits. Quantum deco- 
herence would predict (f> e : <\>^ : (j) T = 1 : 1 : 1, 
independent of source flavor ratios. 

Three-dimensional propagation in structured 
large scale magnetic fields also has considerable 
influence on secondary 7— ray and neutrino fluxes. 
Fig. [2] demonstrates how magnetic fields of /iG 
strength surrounding a UHECR source, for ex- 
ample in a galaxy cluster, can influence the sec- 
ondary GeV-TeV 7— ray fluxes produced by elec- 
tromagnetic cascades initiated by UHECR inter- 
actions with the CMB and infrared background. 
This is the result of simulations with our public 
code CRPropa [35], discussed in Ref. |39l50j . For 
the steep injection spectrum oc E~ 2,7 assumed in 
Fig. [2 the photon flux below a TeV is dominated 
by synchrotron radiation from electron/positron 
pairs produced by protons around the ankle. As a 
consequence, it depends considerably on strength 
and extension of EGMF around the source. 
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